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We analyze various quantum phases of ultracold bosonic atoms in a periodic one dimensional 
optical superlattice. Our studies have been performed using the finite size density matrix renormal- 
ization group (FS-DMRG) method in the framework of the Bose-Hubbard model. Calculations have 
been carried out for a wide range of densities and the energy shifts due to the superlattice potential. 
At commensurate fillings, we find the Mott insulator and the superfluid phases as well as Mott 
insulators induced by the superlattice. At a particular incommensurate density, the system is found 
to be in the superfluid phase coexisting with density oscillations for a certain range of parameters 
of the system. 
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I. INTRODUCTION 

In a seminal paper in 1998, Jaksch et al. [l[ extended 
the work of Fisher et al. @ and predicted the superfluid 
(SF) - Mott insulator (MI) transition in ultracold bosonic 
atoms in an optical lattice. The observation of this tran- 
sition in an optical lattice marked the beginning of the 
experimental studies of quantum phase transitions in ul- 
tracold atoms arising from strong correlations [|[ . It also 
highlighted the possibility for detailed investigations of 
various kinds of strongly correlated quantum systems in 
optical lattices where the different parameters of these 
systems can be exquisitely controlled (3411 ■ 

In this paper, we analyse the various quantum phases 
exhibited by ultracold bosonic atoms in a periodic one 
dimensional optical superlattice using the finite-size den- 
sity matrix renormalisation group (FS-DMRG) method 
for certain parameters of the system. Previously, it 
has been shown using the exact diagonalization, quan- 
tum Monte Carlo and the mean-field decoupling approx- 
imation methods that ultracold bosonic atoms in opti- 
cal super lattices exhibit different phases, with various 
charge-density-wave orders apart from the usual SF and 
MI phases [7H9|]. In the presence of a disordered and 
quasiperiodic potential, it has been shown that the sys- 



tem exhibits apart from SF and MI phases, the quasi- 
Bose glass and the incommensurate charge-density wave 
phases [HI EH- I n our present work, we vary the den- 
sities of the bosonic atoms from commensurate to in- 
commensurate values for a relatively large value of the 
onsite interaction. The purpose of this investigation is to 
predict quantum phases that arise due to the influence 
of the superlattice in the parameter space that we have 
considered. 

The paper is arranged in the following manner. The 
remainder of this section describes the model for bosonic 
atoms in a one dimensional optical superlattice that we 
have chosen in our work. In Section II, the FS-DMRG 
method and the quantities that we have calculated are 
briefly discussed. The results are presented in Section 
III and our conclusions are summarized in Section IV. 

An optical superlattice is formed by the superposition 
of two optical lattices of different frequencies [Fig. (JTJ]. 
Ultracold atoms in an optical superlattice can be de- 
scribed by a modified Bose-Hubbard model where the 
superlattice potential is explicitly included. 
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In the above equation, A; denotes the shift in the energy 
levels for each site due to the superlattice potential. For 
our work, we have considered an optical superlattice with 
the frequency of one optical lattice being double that of 
the other. Hence, the unit cell comprises of two adjacent 
lattice sites as depicted in Fig. (p}. As a result, \ = A; 
V i = odd integers and A.; = 0; V i = even integers. In 
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FIG. 1: (Color online) Superlattice (c) formed by the super- 
position of two optical lattices (a) and (b), frequency of one 
being the double of the other. This results in an energy shift 
of A for alternate sites. 



the first term of Eq. JT]), denotes pair of nearest 

neighbor sites i and j, t denotes the hopping amplitude 
between them with a\(a,i) being the creation (annihila- 
tion) operator which creates (destroys) an atom at site 
i. In the second term of Eq. ([!}, U represents the on-site 
inter-atomic interaction and hi = a|aj is the number op- 
erator. The Hamiltonian is re-scaled in units of hopping 
amplitude t by setting t — 1 so that all the quantities 
becomes dimensionless. 

In the absence of the superlattice potential, the model 
given in Eq. JT]) reduces to the well studied Bose-Hubbard 
model which exhibits a transition from the gapless, com- 
pressible superfiuid (SF) phase to the gapped, incom- 
pressible Mott insulator (MI) phase at integer densi- 
ties However, for non- integer densities, the system 
remains in the superfiuid phase. 

The effect of the superlattice potential is to break the 
translational symmetry of the system. This leads to a 
change in the local (on site) chemical potential. In the 
present case, the local chemical potential changes at al- 
ternate sites in a periodic manner. The SF-MI transi- 
tion has been predicted for such a system at interger 
fillings [9]. Apart from the SF phase, MI phases with 
various fillings at integer and half-integer densities have 
been predicted depending on the relative values of the dif- 
ferent parameters of the system [j| [l(| ■ The MI phase 
is characterized by a fixed number of atoms in a partic- 
ular site or more specifically in a unit cell, and it occurs 
due to the strong onsite interaction between the atoms. 
In the two-period superlattice that we have considered, 
it is possible to have a new class of insulators in which 
alternate lattice sites occupy a fixed number of atoms. 
However unlike the MI phase which is due to strong in- 
teratomic interactions, these insulators arise because of 
the superlattice potential. We therefore refer to these in- 
sulators as superlattice induced Mott insulators (SLMI) . 
If odd sites are occupied by one atom and none in the 
even ones, then we call it as the SLMI-I phase whereas, if 
odd sites are occupied by 2 atoms with none in the even, 



then it is called the SLMI-II phase. 

II. METHOD OF CALCULATION 

To obtain the ground state wave function and the en- 
ergy eigenvalues corresponding to model ([I]) for a sys- 
tem of N bosons in a lattice of length L, with on-site 
interactions and tunnelling as well as the energy shift 
in alternate lattice sites due to the superlattice poten- 
tial, we use the FS-DMRG method with open boundary 
conditions [HI, [l3| . This method is best suited for one- 
dimensional problems and has been successfully used to 
study the Bose-Hubbard model [l3l - [r7j . In our compu- 
tations, we have considered four bosonic states per lattice 
site and the weights of the states neglected in the density 
matrix formed from the left or the right blocks are less 
than 10 -6 [lH. In order to improve the convergence, at 
the end of each DMRG step, we use the finite-size sweep- 
ing procedure given in [lj, EH] • To obtain the ground 
state wave-function \\Pln), and the corresponding energy 
El{N) for densities ranging from 0.24 to 1.25, we start 
with 4 sites and 4 atoms in the FS-DMRG and increase 
both of them by two at each iteration, till we have 24 
atoms. Then we do not increase the number of atoms, 
but increase the number of lattice sites to 100 by adding 
two sites in each DMRG iteration. After each step of 
the iteration, sweeping was done from left to right and 
also right to left across the entire lattice. This process is 
continued till the energy converges. The superlattice po- 
tential breaks the symmetry between the system and the 
environment. Therefore it is necessary to perform the 
calculations on the system and the environment blocks 
separately at each DMRG iteration. Once the length 
L = 100 is reached, we keep it fixed and then increase 
the number of atoms by adding one at a time and per- 
form a complete DMRG sweep for the convergence of 
the energy. This iterative procedure is continued till the 
number of bosons is equal to 125. By this process we 
obtain the wavefunctions and the corresponding energies 
of model([T]) for densities ranging from 0.24 to 1.25. We 
have carried out the calculations for a range of values of 
A starting from to 15, with t and U fixed at 1.0 and 
10.0 respectively. 

Using the ground state wave- function \\Pln), and the 
corresponding energy El(N), the following quantities are 
calculated to identify the various phases exhibited by the 
system. Firstly the chemical potential of a system for the 
density p — N/L is defined as follows: 



The gapped and gapless phases are distinguished from 
the behavior of p as a function of p [18]. To get infor- 
mation about the on-site density distribution in various 
phases exhibited by the system, we calculate the on-site 
density, (n*), defined as, 

{ni) = (W LN \hi\9 LN ). (3) 
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The momentum distribution n(k) and the structure func- 
tion S(k) arc defined as the Fourier transform of the 
single particle density matrix (a p a q ) = l^LN\o/La q \^LN) 
and the density-density correlation function \n p n q ) — 
{^LN\n p n q \^LN), respectively: 

m = ^J2 e%k{p ~ q) H^ ( 4 ) 



S(k) = I J2e ik(p - q) (n P n g ). (5) 
p,q 

The momentum distribution, n(k) will indicate the 
presence of SF phase in the system. The structure func- 
tion will give information about the presence of any den- 
sity wave order in the system, and also the effect of su- 
perlattice potential on the Brillouin zone boundaries. 

To get an idea of the various gapped and gapless phases 
in the system, we plot the density, p vs the average chem- 
ical potential, p. However to draw the phase diagram, 
we need accurate values of the chemical potential, and 
hence we calculate the chemical potential (both p + and 
p~) in the thermodynamic limit L — > oo. [if is defined 
as the difference in the energy of the system of length L 
if one atom is added to the system, whereas p~[ is the en- 
ergy cost when one atom is removed from the system, i.e., 
lit = E L {N+1)-E L (N) and = E L {N)-E L (N-1). 
We plot these chemical potential values versus the recip- 
rocal of the length (1/L), and then extrapolate to length 
tending to infinity (1/L — > 0) to get the values of the 
chemical potential in the thermodynamic limit. We plot 
these chemical potential values for both density half and 
one for different superlattice potential A to get the phase 
diagram. 



SLMFI phase transition as the strength of the superlat- 
tice potential is increased. However, for p = 1 we get 
the MI and the SLMFII phases along with a SF phase 
sandwiched between them. Such a transition at p = 1.0 
had already been predicted in an optical superlattice Q 
using a quantum Monte Carlo approach. However this 
SF phase at p = 1 is not the usual SF phase which arises 
due to a large value of t/U. In fact, in our present case, 
the ratio is very small (t/U ~ 0.1). This SF phase is a 
result of the competition between the superlattice poten- 
tial, A and the on-site repulsive interaction potential, U. 
Moreover, because of the superlattice potential, this SF 
phase will have finite density modulations. 




III. RESULTS AND DISCUSSION 

Before we discuss the details of our results, we first 
summarize the main features of the phase diagram. It is 
well established that bosons in a normal optical lattice, 
described by the Bose- Hubbard model [2[ , exhibit the su- 
pcrfluid to Mott insulator quantum phase transition at 
some critical value of U for integer densities. (For exam- 
ple, for p = 1 this transition occurs at the critical on-site 
interaction U c ~ 3.4 [iHfllj in one-dimension) . For non- 
integer densities, only the superfluid phase exists in the 
ground state. The Mott insulator phase has a finite gap 
and is incompressible. On the other hand, the superfluid 
phase is gapless and compressible. When we include the 
superlattice potential, the energy levels at each site is 
shifted by A^ and this leads to the additional incompress- 
ible, gapped SLMI phase. In model JT]), the SLMI phase 
occurs at commensurate densities, i.e., when p is equal 
to an integer or an half-integer, because the superlattice 
we have considered in our calculations has a periodicity 
of two lattice sites. For p — 1/2, we observe a SF to 



FIG. 2: (Color online) Density p is plotted against the chem- 
ical potential for various values of A at a fixed U = 10 and 
t = 1. 

We first discuss the results for commensurate densities 
and then for incommensurate densities. In the Fig. ([2]) 
we plot p versus p, defined by Eq. ([2]), for a fixed value of 
17 = 10, but A varying from to 15. The gapped phases 
show up in this type of plots as plateaus with the gap Gl 
equal to the width of the plateau, i.e., Gl — — A*7j 
where p, \~ and are the values of the chemical potential 
at the upper and lower knees of the plateau, respectively. 

In Fig. [5] we see for A = 0, there exists only one plateau 
at p = 1.0 correponding to the MI phase, which is ex- 
pected since the value of U is very large compared to 
t [Hj]. This MI phase survives for small values of A. As 
the strength of A increases, two interesting features ap- 
pear in the plots: (i) a new plateau appears at p = 1/2 
and (ii) the width of the plateau at p = 1 decreases. 
This can be clearly seen in the plots of p versus p for 
various A values as shown in Fig. ©. For A = 0.2, the 
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FIG. 3: The chemical potential, p is plotted against density 
p for various values of A, and U — 10. 



gapped phase exists only at p = 1 [Fig. (|3J)] and is gapless 
at p = 1/2. As A is increased beyond 0.8, two gapped 
phases appear, at p = 1/2 and p = 1, respectively. The 
gapped phase at p = 1/2 occurs due to the transition 
from the SF to the SLMI-I phase with one atom per unit 
cell, occupying alternate sites. 

For A values between 0.8 and 9.6, the system exhibits 
two gapped regions, one at p = 1/2 and the other at 
p = 1. The gap at p = 1/2 increases steadily as A in- 
creases and remains finite. On the other hand, for values 
of A > 0.0, the gap at p = 1 decreases continuously and 
ultimately vanishes at around A = 9.6. This kind of be- 
havior for p = 1 is due to the competition between the 
superlattice potential, A and on-site repulsive interaction, 
U. The on-site interaction U tries to impose the MI phase 
in the system whereas A tends to introduce the SLMI-II 
phase. As long as U is greater than A, the MI phase is 
energetically favorable for p — 1 resulting in a finite gap. 
As A becomes comparable to the value of U (here in this 
calculation we have fixed U = 10), neither the MI phase, 
with one boson at every site nor the SLMI-II phase, with 
two bosons at every alternate site is energetically favor- 
able. This leads to vanishing of the energy gap and the 
system becomes a SF. 

When A is larger than the value of U, the superlattice 




FIG. 4: On-site number density plotted against lattice site 
index at density, p = 0.5. 



potential dominates over the on-site interaction. As a 
result, the gap at p = 1 opens up again, which can be 
seen clearly from Fig. ([2]) and ([3]). This gap is due to the 
system being in the SLMI-II phase. The regions other 
than commensurate densities, i.e. p = 1/2 and 1 remain 
in the gapless SF phase for all values of A. 

To confirm the above findings concerning the existence 
of the MI and the SLMI phases, we plot the on-site aver- 
age number density (n.;) with respect to the site index i 
as shown in Fig. (0} and ([5]), for p = 1/2 and 1, respec- 
tively. For values of A less than 0.8, there are no oscilla- 
tions at density equal to half. However, as the value of 
A is increased, the density oscillations begin to set in the 
system, and at higher values of A, a clear {101010 •••} 
occupancy configuration can be seen confirming the pres- 
ence of SLMI-I phase [Fig. • At p = 1, for low values 
of A, the plot [Fig. ([5])] shows a constant number density 
equal to 1.0 implying the MI phase. As A is increased, 
but less than 9.6, the superlattice potential induces some 
oscillations, but the system is still in the MI phase. Even 
at A = 10.0, the oscillations vary from 1.5 (maximum) to 
0.5 (minimum), significantly different from the SLMI-II 
configuration of 12 2 2 ■ ■ ■}. Once A becomes much 




FIG. 5: On-site number density plotted against lattice site 
index at density, p = 1.0. 




larger than U, the SLMI-II configuration is observed. For 
example, at values of A = 14, the number occupancy in 
the sites look very similar to the SLMI-II configuration 
as can be seen from figure [Fig. ©]. 

The momentum distribution is then calculated using 
the Eq.©. The value of n(k) at k = is plotted against 
A for three different values of p (0.5, 0.7 and 1.0). A 
finite value of n(k) at k — is a signature of the SF 
phase in the system. These plots are compared with the 
chemical potential versus p plots, Fig. ([2]), to confirm the 
gapped and gapless phases. In the Fig. ©, for p = 0.5, 
n[k = 0) has a very high value when A is smaller than 
0.8, which signifies the SF phase. As A increases beyond 
0.8, n(k = 0) falls off to values very close to zero, im- 
plying certain phase transition from the SF phase to an 
insulating phase. This insulating phase is the SLMI-I 
as confirmed from the (fii) plots [Fig.©]. For p = 1, 
n(k = 0) starts with values very close to zero [Fig. (J6j)] . 
Due to the very high value of U, the system resides in 
the MI state. As A becomes comparable to U, n(k = 0) 
attains a finite peak value, thus implying the transition 
from the MI to SF phase. This result is in agreement 
with the ones obtained from the chemical potential plots 
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FIG. 7: Structure function, S(k) versus the momentum, k, 
for density 0.5. 



[Figs.© and ©]. As A is increased further, n(k = 0) 
again drops to values very close to zero implying another 
phase transition, but this time the transition being from 
SF to SLMI-II. Hence it can be said that this SF is not 
the usual SF which comes due to t/U ratio, but a con- 
sequence of the competition between the parameter U 
and A. At the incommensurate densities, for example, 
p = 0.7 [Fig©], it is observed that n(k — 0) stays fi- 
nite for all values of A, without any signs of it vanishing, 
showing the presence of SF phase throughout the values 
of A considered in this work. 

Figs. ([7] -[9]) gives the structure function as a function 
of k for various values of density (p = 0.5, 1.0 and 0.7, 
respectively). In the normal optical lattice, the period- 
icity is one lattice site. So for a MI phase, the structure 
function peaks at k = ±2n [19]. However in an optical 
superlattice, the periodicity is doubled. As a result, the 
Brillouin zone is halved. Hence in the MI phase for an op- 
tical superlattice, we expect to observe the peaks of S(k) 
at k = ±7r. It can also be observed that for p = 1, the 
peaks do not appear for lower values of A [Fig. ©]. The 
onset of the peaks actually signify the effects of the super- 
lattice potential, A, in the system Q, and thus changing 
the Brillouin zone boundaries. From Fig. ([7]), we see two 
well defined peaks at k = ±7r for p — 0.5 at a very large 
value of A implying the presence of SLMI-I phase. At 
p = 0.7, we see two peaks at k = ±7r, for A more than 
« 2.0. This is a consequence of some density oscillations 
induced in the system by the superlattice potential. 

The value of structure function at k = 7r is plotted in 
Fig. (|10p for some fixed densities against A to see how 
the effect of superlattice potential evolves in the system. 
For p = 0.5, S(ir) starts from very close to zero, increases 
steadily and then becomes almost constant. For p = 1, 
initially there is no peak in the S(k) plot at k — tt. As A 
is increased, even in the MI phase, we see a finite value of 
the structure function at k = it, indicating the onset of 
the effect of the superlattice potential. Once the A over- 
comes U, the value of S(jr) increases rapidly, showing 
the onset of the SLMI-II phase. It is also noted that at 
A « 10, S(k = 7r) is finite, showing some density modu- 




lations in the system. As seen earlier, at the above value 
of A, the system undergoes a transition from MI to SF 
phase. Thus the SF at p = 1.0 coexists with a finite den- 
sity oscilallation in the system. For p = 0.7, S(ir) is very 
close to zero for low values of A, and then increases as A is 
greater than some value, and then goes on increasing im- 
plying the setting in of some density modulations in the 
system. Hence finite peaks of S(k) at k = n [Fig. ©] for 
p = 0.7 implies the effect of superlattice potential in the 
system bringing about density modulation, although the 
system is in superfluid phase as indicated from n(k = 0) 
plot [Fig. ©] for all values of A. So we see a coexistence 
of SF and some sort of density oscillations in the system 
at the incommensurate density of 0.7. 

The phase diagram in the p - A plane for C/ — 10 is 
presented in Fig. (fTTj). The phase diagram consists of 
SF, MI, SLMI-I and SLMI-II phases depending on the 
density and the superlattice potential. 
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FIG. 10: (Color online) S(k — ir) versus A for three different 
density values, p = 0.5, 0.7 and 1.0. 



IV. CONCLUSION 

We have analyzed the various phases exhibited by a 
system of bosons in an optical superlattice with an unit 
cell consisting of two lattice sites, using the FS-DMRG 
method for a large value of the on-site interaction, U. We 
find that at density, p = 1/2, the system is initially in the 
gapless superfluid phase. But as the value of the super- 
lattice potential, A is increased, a finite gap arises, which 
corresponds to the SLMI-I phase corresponding to the 
configuration [1 1 ...]. At integer density, p = 1, the 
system is initially in the gapped MI phase because of the 
large value of U. But as the superlattice potential value 
is increased, the gap keeps on decreasing, and ultimately 
shrinks to zero when A«f7. This gapless SF phase coex- 
ists with a finite S(k — n) in the system implying density 
oscillations. As the value of A is further increased, the 
gap reopens, corresponding to the SLMI-II phase, with a 
configuration [2 2 0...]. At the incommensurate den- 
sity, p = 0.7, we find the simultaneous existence of the 
SF phase along with density oscillations in the system. 
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